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Abstract. The equation of state for quark matter is derived for a nonlocal, chiral quark model 
within the mean field approximation. We investigate the effects of a variation of the form factors 
of the interaction on the phase diagram of quark matter under the condition of j3- equilibrium and 
^ ' charge neutrality. Special emphasis is on the occurrence of a diquark condensate which signals a 

, phase transition to color superconductivity and its effects on the equation of state. We calculate the 

OO ■ quark star configurations by solving the Tolman- Oppenheimer- Volkoff equations and obtain for the 

I transition from a hot, normal quark matter core of a protoneutron star to a cool diquark condensed 

one a release of binding energy of the order of AMc^ ^ 10^^ erg. We study the consequences of 
antineutrino trapping in hot quark matter for quark star configurations with possible diquark con- 
densation and discuss the claim that this energy could serve as an engine for explosive phenomena. 
A "phase diagram" for rotating compact stars (angular velocity-baryon mass plane) is suggested as 
' a heuristic tool for obtaining constraints on the equation of state of QCD at high densities. It has 

Q-i, a critical line dividing hadronic from quark core stars which is correlated with a local maximum 

^ • of the moment of inertia and can thus be subject to experimental verification by observation of the 

Q ' rotational behavior of accreting compact stars. 
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INTRODUCTION 

Color superconductivity in quark matter UJ is one interesting aspect of recent discus- 
sions devoted to the physics of compact star interiors 10] • Since calculations of the en- 
ergy gap of quark pairing predict a value A ~ 100 MeV and corresponding critical tem- 
peratures for the phase transition to the superconducting state are expected to follow 
the BCS relation Tc = 0.57 A, the question arises whether diquark condensation can 
lead to remarkable effects on the structure and evolution of compact objects. If posi- 
tively answered, color superconductivity of quark matter could provide signatures for 
the detection of a deconfined phase in the interior of compact objects (pulsars, Low- 
mass X-ray binaries) via observations. Hong, Hsu and Sannino have conjectured [3] that 
the release of binding energy due to Cooper pairing of quarks in the course of protoneu- 
tron star evolution could provide an explanation for the unknown source of energy in 
supemovae, hypernovae or gamma-ray bursts, see also |4]. Their estimate of energy re- 
lease did not take into account the change in the gravitational binding energy due to the 
change in the structure of the stars quark core. We reinvestigate the question of a pos- 



sible binding energy release due to a color supercondcutivity transition by taking into 
account changes in the equation of state (EoS) and the configuration of the quark star 
selfconsistently and by including the effects of antineutrino trapping [5]. 

As a first step in this direction we will discuss here the two flavor color supercon- 
ducting (2SC) quark matter phase which occurs at lower baryon densities than the color- 
flavor-locking (CFL) one, see [6, 7]. We will investigate the influence of the formfactor 
of the interaction on the phase diagram and the EoS of dense quark matter under the 
conditions of charge neutrality and isospin asymmetry due to j8 -equilibrium relevant for 
compact stars. 

Finally we consider the question whether the effect of diquark condensation which 
occurs in the earlier stages of the compact star evolution ~ 100 s) [8, 9J at temperatures 
r ~ Tc ~ 20 ^ 50 MeV can be considered as an engine for exposive astrophysical 
phenomena like supernova explosions due to the release of a binding energy of about 
10^^ -^ 10^^ erg, as has been suggested before [|3|,|4j]. 



THERMODYNAMICS OF A NONLOCAL CHIRAL QUARK 

MODEL 

The phase structure of electrically and color neutral quark matter in /3 -equilibrium has 
been studied in [7] and it has been shown that at densities relevant for compact star 
interiors the 2SC phase should be dominant over the CFL phase. The latter one could 
only be stable in the very inner core and thus does not occupy a large enough volume 
in order to cause observable effects. Therefore, we consider in the present work two 
flavor quark matter in with 2SC superconductivity only. The more general case which 
includes the CFL phase does not invalidate the scenario developed in the following and 
will be studied in a subsequent work. We consider the grand canonical thermodynamical 
potential for 2SC quark matter within a nonlocal chiral quark model [10] where in 
the mean field approximation the mass gap and the diquark gap A appear as order 
parameters which can be expressed as in lUlll by 

aq{(l),A;^q,^i,T) = ^^+^-^ dqq^{co[£r{-^q-^i),T] 

+ C0 [£r{llq-lll)j] +C0[er{-llq + lll)j] + 

2 f°° 

+(0[£r{^q + ^l),T]}-^ dqq^{(0[£b{E{q)-[lq)-^I,T] 

+(0 [eb{E{q) + llq) - M/, T]+(0 [£b{E{q) - ^q) + M/, T] 

+(0 [£b{E{q) + pLq) + Ai/, T] } + Q.,ac , (1) 

where we have introduced the quark chemical potential flq = (/i;, + /it/)/2 and the 
chemical potential of the isospin asymmetry /i/ = — /i^) /2 instead of the chemical 
potentials of up and down quark flavors. We neglect here a possible difference between 
the chemical potentials of paired and unpaired colors for the same quark flavor. For 
a more general approach see The factor 2 in the last integral comes from the 



degeneracy of the blue and green colors (£i, = £g). We have introduced the notation 



w[ec,r] = nn[l+exp(-^ 
where the first argument is given by 



+ |. (2) 



e,{x)=x^\+Ai/x\ (3) 

with the color index c = r,g,b and we assume that the green and blue colors are paired 
while the red one remains unpaired, so that we have 

A, = ^(^)A(5,,z, + 5,,g). (4) 

The dispersion relation for unpaired quarks with dynamical mass function m{q) — 
m + g{q)^ is given by 



Ef{q) = sjq^ + m\q) , (5) 

where g{q) denotes the formfactor of the quark interaction, for which we employ the 
following models 

^l(^) = [l + (^/AL)2"]-^ a>l, (6) 
gG{q) = exp(V/AG), (7) 
^njl(^) = 0(1-^/Anjl). (8) 

Depending on the parameter a, the Lorentzian (L) momentum distribution can interpo- 
late between a soft Gaussian (G) formfactor (for a ~ 2) and a hard cutoff (NJL) one (for 
a > 30). The parametrization of the model can be found in Refs. lfl3[floll . The contri- 
bution from the leptons should be added to the quark thermodynamical potential Q.q in 
order to obtain the total one 

le{e,Ve} 

where 

(10) 

is the thermodynamical potential for an ideal gas of massless fermions. The stellar 
matter in the quark core of compact stars consists of u and d quarks, electrons e and 
antineutrinos Vg under the conditions of /3 -equilibrium: d < — > u + + Vg, which in 
terms of chemical potentials reads /ig + /iy^ = — 2jU/ , and charge neutrality: ^riu — ^n^ — 
He — 0, which also could be written as nq + 3nj — 6ne = . The number densities nj, 
j e {u,d,e,Ve,qJ} are defined as 

do. 

^Mj <l>oAo;T 



The conditions for the local extrema of Q.q, correspond to coupled gap equations for the 
two order parameters and A 

do. 

The global minumum of Q.q represents the state of thermodynamical equilibrium from 
which all equations of state can be obtained by derivation. In the following subsections 
we want to comment on aspects of this stellar matter model which turns out to be 
essential for the discussion of quark matter in compact stars. 
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Effect of formfactors 

The nonvanishing of the order parameters or A signals the presence of a phase with 
broken chiral symmetry or color superconductivity, respectively. In Fig. [l] we show the 
resulting phase diagram of quark star matter under the above constraints and neglecting 
the CFL phase which should appear only at such high densities that it will at best occupy 
a negligible volume in the very inner core of a compact star configuration. From Fig.[l] 




FIGURE 1. Phase diagrams for different form factors: Gaussian (solid lines), Lorentzian a = 2 (dashed 
lines) and NJL (dash-dotted). In the upper panel the comparison with the BCS formula for Tc = 0.57 A{T = 
0,jj.g) giflq) is shown for the Gaussian model. 

we see that the softer the formfactor g{q) the lower the critical temperatures and chem- 
ical potentials for the phase transition to quark matter with vanishing order parameters 
or to color superconducting quark matter at low temperatures. It is remarkable that a 
modified BCS relation for the critical 2SC temperature holds ifTol] . 



Antineutrino trapping 



The results for the solution of the gap equations (fT^ are shown in Fig.|3(left panel) for 
different values of the antineutrino chemical potential which is a measure for the density 
of trapped antineutrinos in a hot, young protoneutron star. For values /iy^ > 72 MeV the 
flavor asymmetry becomes large enough to prevent diquark pairing and therefore color 
superconductivity at low densities. Simultaneously, the onset density for quark matter 




FIGURE 2. Left panel: Mass gap (j), diquark gap A and isospin chemical potential /i/ as a function of 
the quark chemical potential jj.^ for different values of the antineutrino chemical potential fiv^ . Solutions 
obey j3 -equilibrium and charge neutrality conditions. Right panel: Pressure vs. energy density for different 
values of the antineutrino chemical potential jiy^. The onset of the appeareance of quark matter is shifted 
to higher energy densities due to antineutrino trapping. 

occurence is shifted to higher densities, see Fig. |2l (right panel). These solutions will 
be used for the discussion of a scenario of hot quark star evolution in the next Section. 
Before that we have to consider the question whether the quark-hadron phase transition 
would occur at too high energy densities so that no quark core could exist in a hybrid 
compact star. 



Quark hadron phase transition 



We construct a quark hadron phase transition at zero temperature using a linear and 
a nonlinear Walecka model for the hadronic phase and perform a Maxwell construction 
for the phase transition. Although it has been claimed that this procedure might be in 
contradiction with the Gibbs conditions for phase equilibrium when more than one 
conserved charge exists in the system |14], a recent investigation including charge 
screening has revealed the opposite |15]. The resulting EoS with a deconfinement 
transition confirms that the lowest critical phase transition (energy) density is obtained 
for the softest (Gaussian) formfactor model, see Fig.|3l 
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FIGURE 3. EoS of compact star matter with quark hadron phase transition. Hadronic EoS: linear (LW) 
and nonUnear (NLW) Walecka model, quark matter: nonlocal separable model with 2SC quark matter and 
different form factors: Gaussian (Gauss), Lorentzian a — 2 (L2) and cutoff (NJL). 

COMPACT STAR CONFIGURATIONS 

What compact star configurations will result from the use of the above EoS? Will a 
stable compact star with quark matter core be obtained? We will answer these questions 
in the present Section by solving the Tolman-Oppenheimer-Volkoff equations for the 
case without rotation and by applying the perturbative method of solution of Einstein 
equations for stationary rigid rotation. First we consider quark stars without hadronic 
shell which can be thought of as the simplified models for a compact star interior and 
after that we discuss hybrid stars. 



Quark stars - engine for explosive phenomena? 

The engine which drives supernova explosions and gamma ray bursts being among the 
most energetic phenomena in the universe remains still puzzling. The phase transiton to 
a quark matter phase may be a mechanism that could release such an amount of energy 
|I16li17.1 . It has been proposed that due to the Cooper instability in dense Fermi gas cold 
dense quark matter shall be in the color superconducting state with a nonvanishing di- 
quark condensate lITsL l^l. The consequences of diquark condensation for the cooling of 
compact stars due to changes in the transport properties and neutrino emissivities have 
been investigated much in detail, see I>19l.20l .8.. .1Q1 . and may even contribute to the ex- 
planation of the relative low temperature of the pulsar in the supernova remmant 3C58 



||2l']. Unlike the case of normal (electronic) superconductors, the pairing energy gap in 
quark matter is of the order of the Fermi energy so that diquark condensation gives con- 
siderable contributions to the equation of state (EoS) of the order of (A/ ji)^. Therefore, 
it has been suggested that there might be scenarios which identify the unknown source of 
the energy of 10^^ erg with a release of binding energy due to Cooper pairing of quarks 
in the core of a cooling protoneutron star [3]. In that work the total diquark condensation 
energy released in a bounce of the core is estimated as (A//i)^Mcore corresponding to 
a few percent of a solar mass, that is 10^^ erg. It has been shown in [10] by solving 
the selfconsistent problem of the star configurations, however, that these effects due to 
the softening of the EoS in the diquark condensation transition lead to an increase in 
the gravitational mass of the star contrary to the naive estimates and that no explosion 
occurs. 

Therefore, we have suggested a new mechanism of energy release fs*] which involves 
a first order phase transition induced by antineutrino untrapping. The trapping of an- 
tineutrinos occurs in hot compact star configurations at temperatures T > \ MeV where 
the mean free path of (anti-)neutrinos becomes smaller than the typical size of a star 
1E2I1 . During the collapse in the hot era of protoneutron star evolution, antineutrinos are 

T = 40 MeV T > 1 MeV ^ _ „ 




)I =72MeV |I =72MeV ^_=0 

\ \ \ 



FIGURE 4. Left two graphs: Quark star cooling by neutrino and photon emission from the surface in 
the case of antineutrino trapping when T > I MeV. Right graph: Antineutrino untrapping and burst-type 
release. 

produced due to the /3 -processes. Since they have a small mean free path, they cannot 
escape and the asymmetry in the system is increased. This entails that the formation of 
the diquark condensate is shifted to higher densities or even inhibited depending on the 
fraction of trapped antineutrinos. As the quark star cools, a two-phase structure will oc- 
cur. Despite of the asymmetry, the interior of the quark star (because of its large density) 
could consist of color superconducting quark matter, whereas in the more dilute outer 
shell, diquark condensation cannot occur and quark matter is in the normal state, opaque 
to antineutrinos for T > 1 MeV. When in the continued cooling process the antineutrino 
mean free path increases above the size of this normal matter shell, an outburst of an- 
tineutrinos occurs and gives rise to an energy release of the order of 10^-^ — 10^"^ erg. This 
untrapping transition is of first order and could lead to an explosive phenomenon. Three 
stages of this scenario of hot quark star evolution are illustrated in Fig. HI In Fig.|5]these 
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FIGURE 5. Three stages of quark star cooling in the phase diagram corresponding to Fig.|3] 

Stages are shown schematically in the quark matter phase diagram. From Fig. El (right 
panel) we can see that the EoS without antineutrinos is softer than with antineutrinos (it 
has a lower pressure at a given energy density) and therefore allows more compact con- 
figurations (Fig.|6l left part of left panel). The presence of antineutrinos tends to increase 
the mass of the star for a given central density (Fig.|6l right part of left panel). 

To estimate the effect of antineutrinos on star configurations we choose a reference 
configuration without antineutrinos with the mass of a typical neutron star Mf = 1.4 Mq 
(see Fig.|6l). The corresponding radius is Rf = 10.29 km and the central density nq = 
10.5 HQ, where = 0. 16 fm^-^ is the saturation density. The configurations with trapped 
antineutrinos and nonvanishing to compare with we choose to have the same total 
baryon number as the reference star: Nb = 1 .48 Nq, where Nq is the total baryon number 
of the sun. For /iy, = 72 MeV we obtain Ma = 1.47 Mq and for /iy, = 150 MeV, 
Mb = 1 .72 Mq. The differences in the radii are Ra—R/ = 0.1 km and RB — Rf = 0.3 km 

and in the central densities — = 0.6 hq and — tiq =2.1 uq, respectively. This is a 
consequence of the hardening of the EoS due to the presence of antineutrinos. The mass 
defect AM,y = Mi — Mf can be interpreted as an energy release if there is a process which 
relates the configurations with Mj and Mf being the initial and final states, respectively. 
In the right panel of Fig. |^ we show limits for the binding energy release as a function 
of the final state mass when conservative values for the antineutrino chemical potential 
during the trapping are chosen. The antineutrino untrapping transition results in a first 
order phase transition which gives rise to an explosive release of the binding energy, as 
required for a scenario which should explain the engine of supernovae or gamma-ray 
bursts i^]. 




FIGURE 6. Left panel: Quark star configurations for different antineutrino chemical potentials /iy^ = 
0, 72, 150 MeV. The total mass M in solar masses Mq is shown as a function of the radius R (left panel) 
and as a function of the central number density ricj in units of the nuclear number density no (right panel). 
Asterisks denote configurations with the same total baryon number. Right panel: Mass defect AM and 
corresponding energy release AE due to antineutrino untrapping as a function of the mass of the final 
state Mf. The shaded region is defined by the estimates for the upper and lower limits of the antineutrino 
chemical potential in the initial state jj,y^ — 150 MeV (dashed-dotted line) and jJ-y^ = 72 MeV (dashed 
hne), respectively. 



Hybid stars 

Once the EoS with a quark-hadron phase transition is defined and the constraints 
of beta equilibrium and charge neutrality are obeyed (see previous Section), the corre- 
sponding hybrid star configurations are obtained from the solution of Einsteins equations 



Is a quark core possible ? 

The answer to this question depends crucially on the employed EoS. Within the setting 
of the present model, stable quark matter cores are possible for a Gaussian formfactor 
model but not for a cut-off one (NJL). This confirms on the one hand conclusions 



previously obtained within other approaches using an NJL model 11241 12511 for quark 
matter but on the other hand presents an alternative quark matter model for which quark 
cores are possible, see Fig. IT] 



Cooling curves and the compact star in 3C58 



In order to obtain the cooling behavior of the presented hybrid star model, we employ 
the program code developed recently so as to describe 2SC quark matter and 
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FIGURE 7. Mass - radius and mass - central energy density curves for hybrid stars. Stable quark matter 
cores are possible for a Gaussian formfactor model but not for a cut-off one (NJL). 

investigate the dependence on the compact star mass ll2lll . The result shown in Fig. [51 
demonstrates that a massive star may have a large enough quark core to entail an enhaced 
cooling behavior in accordance with the recent data point reported by CHANDRA 
measurements of the pulsar in the supernova remnant 3C58. This is, however, not a 
unique feature of a quark matter interior and could also be explained by other exotic 
phases of dense matter Il26ll . 



Population gap for accreting LMXBs 

Our focus is on the elucidation of qualitative features of signals from the high density 
phase transition in the pulsar timing, therefore we use a generic form of an equation of 
state (EoS) with such a transition. We use the polytropic type equation of state for dif- 
ferent values of the incompressibility lilill Ki h{h) = 9 dP/dn at the saturation density, 
see Ref. fr?]. The phase transition between the lower and higher density phases is made 
by the Maxwell construction and compared to a relativistic mean field model consisting 
of a linear Walecka plus dynamical quark model EoS with a Gibbs construction, Fig.|9l 

We introduce a classification of rotating compact stars in the plane of their angular 
frequency Q. and mass (baryon number A^) which we will call phase diagram. In this 
diagram, configurations with high density matter cores are separated from conventional 
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FIGURE 8. Cooling evolution of hybrid stars with different masses. The dots in the inset plot of 
the mass - radius relation correspond to configurations for which the cooling evolution is shown. The 
observational data are taken from 1 26,1 . 



ones by a critical phase transition line. The position and the form of these lines are 
sensitive to changes in the equation of state of stellar matter [5]. 

In Fig.[lO|we display the phase diagrams for the rotating star configurations, which 
correspond to the three model EoS of Fig.|9l These phase diagrams have four regions: 
(i) the region above the maximum frequency Q. > Q.^{N) where no stationary rotating 
configurations are found, (ii) the region of black holes N > N^axi^), and the region of 
stable compact stars which is subdivided by the critical line A^crit(^^) into (iii) the region 
of hybrid stars for > A'crit(^^) where configurations contain a core with a second, high 
density phase and (iv) the region of mono-phase stars without such a core. 

From the comparison of the regional structure of these three different phase diagrams 
in Fig. Uni with the corresponding EoS in Fig. 13 we conclude that there are the following 
correlations between the topology of the lines Nmaxi^) and A/crit(^) and the properties 
of two-phase EoS: 

- The hardness of the high density EoS determines the maximum mass of the star, 
which is given by the line A^niax(^)- Therefore A^max(O) is proportional to the 
parameter KniriH), where nn is the density of the transition to high density phase. 

- The onset of the phase transition line A/crit(0) depends on the density hh and Kt{nL) 
where ni is the density of the transition to the low density phase. 




- The curvature of the lines Nmnx{^) and A'crit(^^) is proportional to the compress- 
ibility of the high and low density phases, respectively. 

Therefore, a verification of the existence of the critical lines A'crit(^) and A^max(^^) by 
observation of the rotational behavior of compact objects would constrain the parameters 
of the EoS for neutron star matter. We have investigated different trajectories of rotating 
compact star evolution in the phase diagram in order to identify scenarios, which result 
in signatures of the deconfinement phase transition lIZTIl . A key evolutionary track 
corresponds to accretion with strong magnetic fields [28]. For this case the Q. first 
decreases as long as the moment of inertia monotonously increases with A^. When 
passing the critical line A^crit(^) for the phase transition, the moment of inertia starts 
decreasing and the internal torque term K^^^^ changes sign. This leads to a narrow dip for 
^{N) in the vicinity of this line. As a result, the phase diagram gets overpopulated for 

N ~ A^crit(^) and depopulated for ~ A^crit(^) up to the second maximum of I{N,Q.) 
close to the black- hole line N^^xi^)- A population gap in the phase diagram of compact 
stars appears as a detectable indicator for hybrid star configurations. 



CONCLUSIONS 

We have investigated the influence of the diquark condensation on the thermodynamics 
of the quark matter under the conditions of /3 -equilibrium and charge neutrality relevant 
for the discussion of compact stars. The EoS has been derived for a nonlocal chiral 
quark model in the mean field approximation, and the influence of different formfactors 




FIGURE 10. Phase diagrams for rotating star configurations corresponding to three model EoS dis- 
cussed in the text. 



(Gaussian, Lorentzian, NJL) has been studied. We have shown that the smoothness of 
the interaction changes the critical temperatures and chemical potentials for the onset of 
the phase transition to lower values. 

The phase transition to color superconducting quark matter from the lower density 
regions at small temperatures (T < 5 ^ 10 MeV) is of first order, while the melting of 
the diquark condensate and the corresponding transition to normal quark matter at high 
temperatures is of second order. The presence of flavor asymmetry due to /3 -equilibrium 
in quark matter does not destroy the diquark condensate since the electron fraction 
Rg/ntotai < 0-01 is too small. The masses of the quark core configurations could be up to 
1.7 Mq and the radii could be up to 1 1 km. 

We have investigated the effects of trapped antineutrinos on the asymmetry and 
diquark condensates in a quark star configurations. By comparing configurations with 
fixed baryon number the release of energy in an antineutrino untrapping transition is 
estimated to be of the order of 10^^ erg. Such a transition is of first order so that 
antineutrinos can be released in a sudden process (burst). This scenario could play an 
important role to solve the problem of the engine of supernova explosions and gamma 
ray bursts. A detailed neutrino transport and cooling calculation should be taken into 
account in a future work. A second antineutrino pulse is suggested as an observable 
characteristics of the present scenario. 

As there are still many unknowns in the picture we have drawn in this contribution 
for the possible effects of quark matter and color superconductivity on compact stars we 
would like to point out that ther will not be one "smoking gun" type signal of quark 
matter but rather the different facets of the picture which emerges for characteristic 



features of compact stars containing quark matter should all match in the big puzzle. 
To bring the pieces together is a task very similar to that of quark gluon plasma search 
in heavy-ion collision experiments. 
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